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ABSTRACT 

A detailed high-resolution spectroscopic analysis is presented for the carbon- 
rich low metallicity Galactic halo object CS 22964-161. We have discovered that 
CS 22964-161 is a double-lined spectroscopic binary, and have derived accurate 
orbital components for the system. From a model atmosphere analysis we show 
that both components are near the metal-poor main-sequence turnoff. Both stars 
are very enriched in carbon and in neutron-capture elements that can be created 
in the s-process, including lead. The primary star also possesses an abundance 
of lithium close to the value of the "Spite-Plateau" . The simplest interpretation 
is that the binary members seen today were the recipients of these anomalous 
abundances from a third star that was losing mass as part of its AGB evolution. 
We compare the observed CS 22964-161 abundance set with nucleosynthesis pre- 
dictions of AGB stars, and discuss issues of envelope stability in the observed 
stars under mass transfer conditions, and consider the dynamical stability of the 
alleged original triple star. Finally, we consider the circumstances that permit 
survival of lithium, whatever its origin, in the spectrum of this extraordinary 
system. 
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1. INTRODUCTION 

The chemical memory of the Galaxy's initial elemental production in short-lived early- 
generation stars survives in the present-day low-mass, low metallicity halo stars showing 
starkly contrasting abundance distributions. Metal-poor stars have been found with order- 
of-magnitude differences in lithium contents, large ranges in a-element abundances (from 
[Mg,Si,Ca,Ti/Fej^ < to ~+l), non-solar Fe-peak ratios, and huge bulk variations in 
neutron-capture (n-capture) abundances. 

One anomaly can be easily spotted in medium-resolution (R = A/AA ^ 2000) spec- 
troscopic surveys of metal-poor stars: large star-to-star variations in CH G-band strength, 
leading to similarly large carbon abundance ranges. Carbon-enhanced metal-poor stars 
(hereafter CEMP = [C/Fe] > +1) are plentiful at metallicities [Fe/H] < -2, with esti- 
mates of their numbers ranging from ~14% (Cohen et al. 2005) to ^21% (Lucatello et al. 
2006). The large carbon overabundances are mainly (but not always) accompanied by large 
n-capture overabundances, which are usually detectable only at higher spectral resolution 
(R > 10,000). In all but one of the CEMP stars discovered to date, the n-capture abun- 
dance pattern has its origin in slow neutron-capture synthesis (the s-process). The known 
exception is CS 22892-052 {e.g., Sneden et al. 2003a and references therein), which has 
[C/Fe] ~ -|-1 but an n-capture overabundance distribution clearly consistent with a rapid 
neutron-capture (the r-process) origin. Properties of CEMP stars have been summarized 
recently in a large-sample high-resolution survey by Aoki et al. (2007). They conclude in 
part that the abundances of carbon and (n-capture) barium are positively correlated (their 
Figure 6), pointing to a common nucleosynthetic origin of these elements in many CEMP 
stars. 



^This paper includes data gathered with the 6.5ni MageUan and 2.5m du Pont Telescopes located at Las 
Campanas Observatory, Chile. 

^ We adopt the standard spectroscopic notation (Heifer, Wallerstein, & Greenstein 1959) that for elements 
A and B, log e(A) = logio(NA/NH) + 12.0, and [A/B] = logio(NA/NB)* - logiolNA/Ne)©. Also, metallicity 
is defined as the stellar [Fe/H] value. 
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CS 22964-161 was first noted in the "HK" objective prism survey of low metallicity halo 
stars (Beers, Preston, & Shectman 1992, hereafter BPS92). Using as metallicity calibra- 
tion the strength of the Ca II K line, BPS92 estimated [Fe/H] ~ —2.62. They also found 
CS 22964-161 to be one of a small group of stars with unusually strong CH G-bands (see their 
Table 8). Recently, Rossi et al. (2005) analyzed a moderate resolution spectrum of this star. 
From the extant BVJK photometry and distance estimate they suggested that CS 22964-161 
is a subgiant: Teg = 5750 K and log g = 3.3. Two metallicity estimates from the spectrum 
were in agreement at [Fe/H] = —2.5, and three independent approaches to analysis of the 
overall CH absorption strength suggested a large carbon abundance, [C/Fe] = +1.1. 

We observed CS 22964-161 as part of a high resolution survey of candidate low- metallicity 
stars selected from BPS92. When we discovered that the star shows very strong features of 
CH and n-capture species Sr II and Ba II, similar to many binary blue metal-poor (BMP) 
stars (Sneden, Preston, & Cowan 2003b), we added it to a radial velocity monitoring pro- 
gram. Visual inspection of the next observation showed two sets of spectral lines, and so an 
intensive monitoring program was initiated on the du Pont and Magellan Clay telescopes at 
Las Campanas Observatory. In this paper we present our orbital and abundance analyses of 
CS 22964-161. Radial velocity data and the orbital solution are given in ^ and broadband 
photometric information in ^ We discuss the raw equivalent width measurements for the 
combined-light spectra and the extraction of individual values for the CS 22964-161 primary 
and secondary stars in ^ Determination of stellar atmospheric parameters is presented in 
^ followed by abundance analyses of the individual stars in ^ and of the CS 22964-161 
system in §71 Interpretation of the large lithium, carbon, and s-process abundances in pri- 
mary and secondary stars is discussed in ^ Finally, we speculate on the nature of former 
asymptotic giant branch (AGB) star that we suppose was responsible for creation of the 
unique abundance mix of the CS 22964-161 binary in fj9l 



2. RADIAL VELOCITY OBSERVATIONS 

We obtained spectroscopic observations of CS 22964-161 with the Clay 6.5-m MIKE 
(Bernstein et al. 2003) and du Pont 2.5-m echelle (R ~ 25,000) spectrographs. Properties 
of the two spectrographs are presented at the LCO websit^. The Clay MIKE data have 
R = X/6X and continuous spectral coverage for 3500 A < A < 7200 A. The du Pont data 
have R ^ 25,000, and out of the large spectral coverage of those data we used the region 
4300 A < A < 4600 Afor velocity measurements. Exposure times ranged from 1245 to 3500 



•^http:/ /www.lco.cl 
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seconds on the Clay telescope and 3000 to 4165 seconds on the du Pont telescope. The 
observations generally consist of two exposures flanked by observations of a thorium-argon 
hollow-cathode lamp. 

The Magellan observations were reduced with pipeline software written by Dan Kelson 
following the approach of Kelson (2003, 2006). Post-extraction processing of the spectra 
was done within the IRAF ECHELLE packaged The du Pont observations were reduced 
completely with IRAF ECHELLE software. 

Velocities were initially measured with the IRAF FXCOR package using MIKE obser- 
vations of HD 193901 as a template, and a preliminary orbit was derived. The three MIKE 
observations of CS 22964-161 obtained at zero velocity crossing (phase ^ 0.11) were then 
averaged together to deflne a new template, hereafter called the syzygy spectrum. This spec- 
trum has a total integration time of 10215 sec and S/N > 160 at 4260 A. We remeasured the 
velocities using this new template with the TODCOR algorithm (Zucker & Mazeh 1994). 
The cross-correlations covered the wavelength interval 4130 A < A < 4300 A. The syzygy 
spectrum was also used extensively in our abundance analysis (^HD- 

The radial velocity data were flt with a non-linear least squares solution. We chose 
to fit to only the higher resolution and higher S/N Magellan data, using the du Pont data 
to confirm the orbital solution. The observations are presented in Table [1] which lists the 
heliocentric Julian Date (HJD) of mid-exposure, the velocities of the primary and secondary 
components, and the orbital phases of the observations. The adopted orbital elements are 
listed in Table [2] and the adopted orbit is plotted in Figure [H Of particular importance to 
later discussion in this paper are the derived masses: Mp sin^ i = 0.773 ± 0.009 Mq and 
Ms sin^ i = 0.680 ± 0.007 Mq; the orbital inclination cannot be derived from our data. 

We will discuss the velocity residuals to the orbital solution further in §21 Hereafter, 
references to individual spectra will be by the HJD of Tabled! e.g., the spectrum obtained 
with Magellan MIKE on HJD 13817.90625 will be called "observation 13817". 

3. PHOTOMETRIC OBSERVATIONS 

Preston, Shectman, & Beers (1991) list V = 14.41, B-V = 0.488, andU-B = -0.171 
for CS 22964-161 based on a single photoelectric observation on the du Pont telescope. New 
CCD observations of this star were obtained with the du Pont telescope on UT 28 May 



''IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under cooperative agreement with the NSF. 
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2006. The data were calibrated with observations of the Landolt standard field Markarian 
A (Landolt 1992). Observations of the standards and program object were taken at an 
airmass of ~1.12, and we used standard extinction coefficients in the reductions. We derived 

V = 14.43 and B — V = 0.498 for CS 22964-161. Old and new photometric data are consistent 
with the typical ~1% observational uncertainties in magnitudes, so we adopt final values of 

V = 14.42 and B -V = 0.49. 

The ephemeris given in Table [2] predicts a primary eclipse at HJD2454315.045 (2 Au- 
gust 2007, UT 13.08 hours) and a secondary echpse at HJD 2454344.453 (31 August 2007, 
UT 22.87 hours). CS 22964-161 was monitored on the nights of 2/3 August 2007 and 31 
August /I September using the CCD camera on the Swope telescope at Las Campanas. Ob- 
servations were obtained approximately every hour with a V filter. No variations in the V 
magnitude of CS 22964-161 in excess of 0.02 magnitudes were detected. 



A double-line binary spectrum is a complex, time-variable mix of two sets of absorption 



observed equivalent widths EWo to be those measured in the combined-light spectra, and true 
equivalent widths EWt to be those of each star in the absence of its companion's contribution. 
If the primary-secondary velocity separation is large, one can attempt to measure observed 
EWs of primary and secondary stars independently. If this can be accomplished then true 
EWs can be computed from knowledge of the relative fiux levels of the stars. In practice 
the derivation of true EWs is complex and subject to large uncertainties. 

In Figure [2] we illustrate the difficulties in deducing observed EWs for each star from 
the observed CS 22964-161 composite spectra. Recall that in the syzygy spectrum (the co- 
addition of three individual observations), there is no radial velocity separation, or AVr ~ km s 
The primary and secondary spectral lines coincide in wavelength, producing the simple spec- 
trum shown at the bottom of Figure [2J In contrast, the middle spectrum in this figure was 
generated from spectra in which the velocity separation was large enough for primary and 
secondary lines to be resolved (hereafter called velocity-split spectra). We co-added three of 
these spectra that have similarly large velocity differences, AVr ~ 60 km s^^ (observations 
13580, 60.12 km s'^ 13817, 58.09 km s'^; and 13832, 62.40 km s"^ Table The weaker, 
redshifted secondary absorption lines are obvious from comparison of this and the syzygy 



4. EQUIVALENT WIDTH DETERMINATIONS 



lines and two usually unequal continuum flux levels. Following Preston 




A similar kind of analysis of the very metal-poor binary CS 22876-032 has been discussed by Norris, 
Beers, & Ryan (2000). 
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spectrum. 

Identification of tlie secondary absorption lines is clarified by inspection of tfie top 
spectrum in Figure [H We created tliis artificial spectrum to mimic tfie appearance of tfie 
secondary by diluting tfie observed syzygy spectrum witli addition of a constant tliree times 
larger than the observed continuum, and shifting that spectrum redward by 60 km s~^. In 
this line-rich wavelength domain, the velocity shift yields a few cleanly separated primary 
and secondary absorption features. For example, 4199.11 A Fe I in the primary becomes 
4199.95 A in the secondary, and observed EWs of both stars can be measured. More 
often however, secondary lines are shifted to wavelengths very close to other primary lines, 
destroying the utility of both primary and secondary features. An example of this is the 
Fe I 4198.33 A line, which at AVr ~ 60 km s~^ becomes 4199.17 A for the secondary star, 
which contaminates the Fe I 4199.11 A line of the primary star. The various blending issues, 
combined with the intrinsic weakness of the secondary spectrum, yield few spectral features 
with clean EW values for both primary and secondary stars. 

4.1. Equivalent Widths from Comparison of Syzygy and Velocity-Split Spectra 

In view of the blending issues outlined above, we derived the observed EWs through 
a spectrum difference technique. On the mean syzygy spectrum, we measured EWo^tot = 
EWo,p + EWo,s, where subscripts o denotes an observed EW, and p and s denote primary 
and secondary stars. On five spectra with large primary- secondary velocity separations 
(13580, 13585, 13587, 13817, 13834) we measured the primary star lines, EWo^p. These five 
independent values were then averaged, and then the secondary star's values were computed 
as EWo^s = EWo^tot— < EWo^p >■ In this procedure we attempted to avoid primary hues 
that would have significant contamination by secondary lines in the velocity-split spectra. 

The observed EW values are given in Table [3l along with the line excitation potentials 
and transition probabilities. In this table we also include the atomic data for lines from which 
abundances ultimately derived from synthetic spectrum rather than EW computations. For 
one estimate of the uncertainty in our EW measurement procedure, we computed standard 
deviations for the primary-star measurements of each line aEWo,p- The mean and median of 
these values for the whole line data set were 5.2 and 4.1 mA, respectively. The uncertainties 
in EWo^tot are expected to be smaller because the syzygy spectrum is the mean of three 
individual observations and is thus of higher S/N. Therefore we take 4—5 mA as an estimate 
of the uncertainty in EWo,s values determined in the subtraction procedure. 

For a few very strong transitions the secondary lines are deep enough to cleanly detect 
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when they split away from the primary hnes. We have employed these to assess the reliability 
of the EWs derived by the subtraction method described above. In Figure [3] we show a 
comparison of EWo^s values given in Table [3] for seven lines that we measured on up to six 
velocity-split spectra (the five named in the previous paragraph plus 13832, which was not 
used in the subtraction procedure). Taking the means of the EWo,s measurements of each 
line and comparing them to the subtraction-based values of Table [3], the average difference for 
the seven lines is 0.8 mA with a scatter a = 3.5 mA (the median difference is 0.6 mA). This 
argues that in general the individual EWo^s values agree with the EWo^s subtraction-based 
ones. 

Derivation of true EWs depends on knowledge of the relative luminosities of primary 
and secondary stars. Formally, from Preston's (1994) equations 4—5, we have EWp^t = 
EWp^oi^ + ^s/lp), and EWs^t = -^^5,0(1 + Ip/h), where subscript t represents the true EW , 
and / denotes an apparent luminosity. The luminosity ratios will be wavelength-dependent if 
the two stars are not identical in temperature. For the entire line data set, ignoring weaker 
primary lines (those with EWp^o < 25.0 mA), we calculated a median equivalent width ratio 
{EWp^o/ EWs^o) = 5.2. Inspection of the velocity-split spectra suggested that relative line 
strengths in the secondary spectrum were not radically different from those of the primary 
spectrum. Therefore we concluded that the spectral types of the stars are not too dissimilar 
and therefore we concluded that to first approximation EWp^t ~ EWg^t- This assumption 
then leads to EWp^o/ EWs,o ~ ip/h ~ 5. 

To account roughly for the small derived temperature difference (see f|5]) we finally 
adopted Ip/ls = 5.0 in the photometric V bandpass (A ~ 5500 A), and increased the ratio 
linearly by a small amount with decreasing wavelength. Thus in the B bandpass (A ~ 4400 A) 
we used Ip/ls = 5.6. Final true EW values using this prescription are given in Table [31 The 
correction factors between observed and true EWs were approximately 1.2 and and 6.0 for 
primary and secondary stars in the V spectral region, and the disparity in these factors is 
larger at B. Clearly the EWp^t values of Table [3] are much more reliably determined than 
the EWs^t ones. 



5. STELLAR ATMOSPHERE PARAMETERS 

5.1. Derivation of Parameters 

We used the observed B — V color and minimum masses for the binary together with 
the Victoria- Regina stellar models (VandenBerg, Bergbusch, & Dowler 2006), to estimate 
initial model atmospheric parameters for the component stars of CS 22964-161. We adopted 
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B — V = 0.49 from ^and E{B — V) = 0.07 (BPS92; a value we also estimate employing the 
dust maps of Schlegel, Finkbeiner, & Davis 1998) to obtain {B — V)o = 0.42. We assumed 
initially that sin^ i ~ 1.0 for the binary orbit and thus Mp = 0.773 Mq and Mg = 0.680 Mq, 
as derived in §2J We interpolated the Victoria-Regina models computed for [Fe/H] = —2.31, 
[a/Fe] = +0.3, and Y = 0.24, to obtain evolutionary tracks for these masses. These tracks 
were used to derive B — V colors and luminosity ratios for the component stars. We adopted 
starting values of Tgfj, log g, and Ip/lg where the tracks give B — V = 0.42 for the combined 
system. These values were Teg p = 6050 K, T^fi^s = 5950 K and log gp = 3.6, log gs = 4.2. 
The luminosity ratios using these parameters were Ip/lg (B) = 8.65 and Ip/ls {V) = 8.39, 
somewhat larger than implied by our spectra. 

Final model atmospheric parameters were found iteratively from the EW data for the 
two stars. We employed the LTE line analysis code MOOG@ (Sneden 1973) and interpolated 
model atmospheres from the Kurucz (1998)0 grid computed with no convective overshoot 
(as recommended by Castelli, Gratton, & Kurucz 1997, and by Peterson, Dorman, & Rood 
2001). 

We began by using standard criteria to estimate the model parameters of the CS 22964-161 
primary: (a) for Tcs, no trend of derived Fe I individual line abundances with excitation po- 
tential; (b) for ft, no trend of Fe I abundances with EW; (c) for log g, equality of mean Fe I 
and Fe II abundances (for no other element could we reliably measure lines of the neutral 
and ionized species); and (d) for model metallicity [M/H], a value roughly compatible with 
the Fe and a abundances. These criteria could be assessed reliably for the primary star 
because it dominates the combined light of the two stars. In the top panel of Figure H] we 
illustrate the line-to-line scatter and (lack of) trend with wavelength of the primary's Fe I 
and Fe II transitions. With iteration among the parameters we derived {Tes^p, log gp, ftp, 
[M/H]p) = (6050±100 K, 3.7±0.2, 1.2±0.3 km s-\ -2.2±0.2). 

To estimate model parameters for the secondary we assumed that derived [Fe/H] metal- 
licities and abundance ratios of the lighter elements (Z < 30) in the primary and secondary 
stars should be essentially identical if they were formed from the same interstellar cloud. 
Iteration among several sets of {Tes,s, log gs) pairs was done until [Fe/H]^ ~ [Fe/H]p and Fe 
ionization equilibrium for the secondary was achieved. Given the weakness of the secondary 
spectrum and the resulting large correction factors used to calculate EWs,t from EWs,o val- 
ues, it is not surprising that the spectroscopic constraints on the secondary parameters were 
weak. This is apparent from the adopted Fe abundances displayed in the bottom panel of 



Available at: http://verdi.as.utexas.edu/moog.html 



^ Available at: http://kiu'ucz. harvard.edu/ 
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Figure HI The o values for individual line abundances were about three times larger for the 
secondary than the primary (Table Hj), and number of transitions was substantially smaller 
(e.g'., we measured seven Fe II lines for the primary but only three for the secondary). 

We adopted final model parameters of (Teg^s, log fts, [M/H]s) = (5850 K, 4.1, 
0.9 km s~^, —2.2). Model uncertainties for the secondary were not easy to estimate and, 
of course, were tied to our opening assumption that the two stars have identical overall 
metallicities and abundance ratios. Thus with [M/H]p = [M/H]s, the uncertainties in Teg, 
log (7, and ft of the secondary are approximately double their values quoted above for the 
primary. Thus caution is obviously warranted in interpretation of the model parameters of 
the CS 22964-161 secondary. 

5.2. Comparison to Evolutionary Tracks 

The well-determined Tefj p and log values can be used with the mass and luminosity 
ratios of the stars to provide an independent estimate of Teg^s and log gs- Standard relations 
L oc R^T^s and g oc MjW lead to 

log(L,/Lp) = log(/,//p) = log(M,/Mp) + 4 X log(Teff,,/Teff,p) - Xog^g^l g^) 

Taking an approximate average luminosity ratio to be l^jls — 5.3, adopting Mg/Mp = 1.15 
(Table [2]) and assuming Tcs,p = 6050 K and log gp = 3.7 from above, we get a predicted 
temperature-gravity relationship log gs — 4xlog(Teff,s) — 10.8 . 

In Figure [5] we show the very metal-poor main sequence turnoff region of the H-R 
diagram in (log T^g, log g) units. The Victoria- Regina evolutionary tracks (VandenBerg 
et al. 2006; [Fe/H] = -2.31, [a/Fe] = +0.3, and Y = 0.24) discussed in are again 
employed. However, the unknown binary orbital inclination of CS 22964-161 cannot be 
ignored here. Therefore we have plotted the tracks for three pairs of masses corresponding 
to assumed sin i values of 90°, 80°, and 75°. Also plotted is a straight line representing the 
CS 22964-161 secondary star temperature-gravity equation derived above. 

The CS 22964-161 primary and secondary and (Tes, log g) positions are indicated with 
filled circles in Figure [51 We also add data indicated with open triangles from CEMP 
high resolution spectroscopic studies for C-rich stars of similar metallicity, taken here to 
be [Fe/H] = -2.4 ± 0.4. The studies include those of Aoki et al. (2002), Sneden et al. 
(2003b), and Cohen et al. (2006). If Tgg = 5850 K from the spectroscopic analysis then 
the temperature-gravity relationship from above predicts log gs = 4.3 (indicated by a filled 
square in the figure). Our spectroscopic value of log gs = 4.1 lies well within the uncertainties 
of both estimates. 
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Given the apparently anomalous position of the CS 22964-161 secondary in Figure O it 
is worth repeating its abundance analysis using atmospheric parameters forced to approx- 
imately conform with the evolutionary tracks. This is equivalent to attempting a model 
near the high temperature, high gravity end of its predicted relationship shown in the fig- 
ure. Therefore we computed abundances for a model with parameter set (Tes,s, log gs, 
[M/H]s) = (6300 K, 4.5, —2.2). Assumption of a secondary microturbulent velocity Vts = 
Vtp = —1.2 km s~^, yields log e(Fe) = 5.45 and 5.27 from Fe I and Fe II lines, respectively. 
These values are substantially larger than the mean abundances for primary and secondary 
given in Table HI log e(Fe) = 5.11. This is in agreement with expectations of a larger derived 
metallicity from the 350 K increase in T^s^s for this test. However, the Fe I line abundances 
for this hotter model also exhibit an obvious trend with EW. This problem could be cor- 
rected by increasing Vt to 2.4 km s~^, and then we get log e(Fe) = 5.18 and 5.15 from Fe I 
and Fe II lines, very close to our final adopted Fe abundances for the secondary. However, 
it is difficult to reconcile such a large Vt with the much smaller value determined with more 
confidence in the primary, as well as standard values determined in many literature studies 
of near-turnoff stars. 



6. ABUNDANCES OF THE INDIVIDUAL STARS 

With the EW data of Table [3] and the interpolated model atmospheres described in 
^ we determined abundances of a few key elements whose absorption lines are detectable 
in both primary and secondary observed spectra. These abundances are given in Table |H 
They suggest that in general the abundance ratios of all elements in the two stars agree 
to within the stated uncertainties. Both stars are relatively enriched in the a elements: 
<[Mg,Ca,Ti/Fe]p > ~ +0.4 and <[Mg,Ca,Ti/Fe]^ > ^ +0.3. Both stars have solar-system 
Ni/Fe ratios, probably no substantial depletions or enhancements of Na (special comment 
on this element will be given in and large deficiencies of Al. All these abundance ratios 
are consistent with expectations for normal metal-poor Population II stars. 

More importantly, we find very large relative abundances of C, Sr, and Ba ([X/Fe] = +0.5 
to +1.5; Table H]) in both primary and secondary stars of the CS 22964-161 system. For Sr 
and Ba abundances we first used the EW subtraction technique, which suggested roughly 
equal abundances of these elements in both stars. We confirmed and strengthened this 
result through synthetic spectrum computations of the strong Sr II 4077.71, 4215.52 A and 
the Ba II 4554.04 A lines in the six velocity-split spectra. To produce the binary syntheses 
we modified the MOOG line analysis code to compute individual spectra for primary and 
secondary stars, then to add them after (a) shifting the secondary spectrum in wavelength 
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to account for the velocity difference between the stars, and (b) weighting the primary and 
secondary spectra by the appropriate luminosity ratio. 

In Figure [HI we show the resulting observed/synthetic binary spectrum match for the 
Sr II 4077.71 A line in observation 13817. The absorption spectrum of the secondary star 
is shifted by +58.1 km s~^ (+0.79 A), in agreement with the observed feature in Figure [61 
The depth of the Sr II line in the secondary is weak, as expected due to the luminosity 
difference between the two stars {Ip/lg — 6) at this wavelength. Note also in Figure the 
relative insensitivity of the feature to abundance changes for the secondary, even with the 
large (±0.5 dex) excursions in its assumed Sr content. This occurs because the Sr II line 
is as saturated in the secondary as it is in the primary, but the central intensity of the 
unsmoothed spectrum is ~0.2 of the continuum. Thus after the secondary's spectrum is 
diluted by the much larger flux of the primary, a weak line that is relatively insensitive to 
abundance changes ensues in the combined spectrum (and a naturally weaker line in the 
secondary simply becomes undetectable in the sum). 

In Figure [7] we illustrate the appearance of a small portion of the CH G-band A^A — 
X'^U'^ Q-branch bandhead in the observed combined-light spectra with a large velocity split. 
This portion of the G-band has a sharp blue edge; the central wavelength of the first line is 
4323.0 A. With AVr = +58.1 km s~^, the left edge of the secondary's bandhead begins at 
4323.8 A, as indicated in the bottom panel of Figure [71 In the top panel we show attempts 
to match the observed bandhead with synthetic spectra that include only CH lines of the 
primary star. When the spectral interval A < 4324.8 A is fit well, observed absorption is 
clearly missing at longer wavelengths in the synthetic spectrum. This is completely solved 
by the addition of the secondary's CH bandhead, at a comparable C abundance level to that 
of the primary, as shown in the bottom panel. 

It is very difficult to derive reliable abundances in the CS 22964-161 secondary even for 
the strong features illustrated here. Nevertheless, it is clear that both components of this 
binary have substantial overabundances of C and n-capture elements Sr and Ba. Within the 
uncertainties of our analysis, the overabundance factors for these elements appear to be the 
same. Enhanced C accompanies s-process synthesis of n-capture elements during partial He- 
burning episodes of low/intermediate-mass stars, and the joint production of these elements is 
evident in the observed abundances of a number of BMP stars such as CS 29497-030 (Sneden 
et al. 2003b; Ivans et al. 2005). However, a substantial C overabundance has also been seen 
in the r-pro cess-rich star CS 22892-052 (Sneden et al. 2003a). The Sr and Ba abundances 
determined to this point cannot distinguish between possible n-capture mechanisms that 
created the very heavy elements in CS 22964-161. 
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7. ABUNDANCES FROM THE SYZYGY SPECTRUM 

To gain further insight into the n-capture element abundance distribution we returned 
to the higher S/N mean syzygy spectrum. Preston et al. (2006b) argued that the relative 
strengths of La II and Eu II lines can easily distinguish s-process dominance (stronger La 
features) from r-process dominance (stronger Eu); see their Figure 1. In panels (a) and (b) of 
Figure [8] we show the same La and Eu lines discussed by Preston et al.; the greater strength 
of the La feature is apparent. 

We computed a mean EW from La II hues at 3988.5, 3995.7, 4086.7, and 4123.2 A, and 
a mean EW from Eu II lines at 3907.1, 4129.7, and 4205.1 A for a few warm metal-poor stars 
with n-capture overabundances. The resulting ratio < -EW^La >/< EWeu > — 0.5 for the 
r-pro cess-rich red horizontal-branch star CS 22886-043 (Preston et al. 2006a), ^ 2.8 for the 
s-process-rich RR Lyrae TY Gruis (Preston et al. 2006a), and ~ 2.7 for the r+s BMP star 
CS 29497-030 (Ivans et al. 2005). For CS 22964-161 we found < EWia >/< EWeu > ^ 2.7, 
an unmistakable signature of an s-process abundance distribution. 

For a more detailed n-capture element distribution for CS 22964-161 we computed syn- 
thetic spectra of many transitions in the syzygy spectrum. We used the same binary syn- 
thesis version of the MOOG code that was employed for the velocity-split spectra shifted 
to AVr = km s~^ (as described in For these computations it was also necessary to 
assume that log e(X)p = log e{X)s for all elements X. We also derived abundances of a few 
lighter elements of interest with this technique. It should be noted that since the primary 
star is 4-6 times brighter than the secondary, abundances derived in this manner mostly 
apply to the primary. 

In Table [5] we give the abundances for the CS 22964-161 binary system. For elements 
with abundances determined for the individual stars as discussed in ^ we also give estimates 
of their "system" abundances in this table. These mean abundances were computed from 
the entries in Table HI giving both the abundances and their uncertainties (a) of the primary 
star five times more weight than those of the secondary. For elements with abundances 
derived from synthetic spectra of the syzygy spectrum, the abundances and a values are 
means of the results from individual lines, wherever possible. For several of these elements 
only one transition was employed. In these cases we adopted a = 0.20 or 0.25, depending 
on the difficulties attendant in the synthetic/observed spectrum matches. In the next few 
paragraphs we discuss the analyses of a few species that deserve special comment. 

Li I: The resonance transition at 6707.8 A was easily detected in all CS 22964-161 
spectra, with EWo^tot — 24.5 mA from the syzygy data. Our synthesis included only ^Li, 
but with its full hyperfine components. The resolution and S/N combination of our spectra 
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precluded any meaningful search for the presence of ^Li. Reyniers et al. (2002) have shown 
that the presence of a Ce II transition at 6708.09 A can substantially contaminate the Li 
feature in n-capture-rich stars. However, in our spectra the Ce line wavelength is too far 
from the observed feature, and our syntheses indicated that the Ce abundance would need 
to be about two orders of magnitude larger than our derived value (Table [5]) to produce 
measurable absorption in the CS 22964-161 spectrum. 

We attempted to detect the secondary Li I line in two different ways. First, we applied 
the subtraction technique (Q to this feature. However, the six velocity-split spectra yielded 
EWo,p = 23.5 mA with a = 3.5 mA (consistent with typical uncertainties for this technique). 
Thus the implied EWq^s = 1.0 mA is consistent with no detection of the secondary's line. 
Second, we co-added the velocity-split spectra after shifting them to the rest system of the 
secondary. A very weak {EWo,s — 4 mA) line appears at A = 6707.9 A, about 0.1 A redward 
of the Li I feature centroid. We computed synthetic spectra assuming that the observed 
line is Li, and derived log e(Li) ~ +2.0 with an estimated uncertainty of ±0.2 from the 
observed/synthetic fit. This abundance is fortuitously close to the system value, but we 
do not believe that the present data warrant a claim of Li detection in the CS 22964-161 
secondary. 

Na I: The 5682, 5688 A doublet, used in abundance analyses when the D-lines at 5890, 
5896 A are too strong, is undetectably weak in CS 22964-161. A synthetic spectrum match 
to the syzygy spectrum suggests that log e(Na) < 3.9. The D lines are strong, ~40% deep. 
Unfortunately, their profiles appear to be contaminated by telluric emission. If we assume 
that the line centers and red profiles are unblended then log e(Na) ~ 3.9. 

Given the great strengths of the D-lines in the syzygy spectrum (dominated by the 
primary star) we returned to the velocity-split spectra and searched for the secondary Na D 
lines. Absorptions a few percent deep at the expected wavelengths were seen in all spectra. 
Co-addition of these spectra after shifting to the rest velocity system of secondary yielded 
detection of both D-lines, and we estimated log e(Na) ~ 3.8. A problem with our analyses 
of the D-lines in all spectra was the lack of telluric H2O line cancellation; we did not observe 
suitable hot, rapidly rotating stars for this purpose. However, any unaccounted-for telluric 
features would drive the derived Na abundances to larger values. Thus we feel confident that 
Na is not overabundant in CS 22964-161. 

CH: Carbon was determined from CH G-band features in the wavelength range 4260- 
4330 A, with a large number of individual lines contributing to the average abundance. Our 
CH line list was taken from the Kurucz (1998) compendium. The solar C abundance listed 
in Table [5] was determined with the same line list (see Sneden et al. 2003b), rather than 
adopted from the recent revision of its abundance from other spectral features by AUende 
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Prieto, Lambert, & Asplund (2002). 

Si I: The only detectable line, at 3905.53 A, suffers potentially large CH contamination, 
as pointed out by, e.g., Cohen et al. (2004). These CH lines are very weak in ordinary warm 
metal-poor stars (Preston et al. 2006a), but are strong in C-enhanced CS 22964-161. We 
took full account of the CH in our synthesis of the Si I line. Note that the secondary's Si I 
line can be detected in the velocity-split spectra, but it is always too blended with primary 
CH lines to permit a useful primary /secondary Si abundance analysis. 

Species with substructure: Lines of Sc II, Mn I, Y II, and La II have significant hy- 
perfine subcomponents, which were explicitly accounted for in our syntheses (each of these 
elements have only one naturally-occurring isotope). For Ba II and Yb II (whose single 
spectral feature at 3694.2 A is shown in panel (c) of Figure [HI both hyperfine and isotopic 
substructure were included. We assumed an s-process distribution of their isotopic fractions: 
/(i34Ba) = 0.038, fC^^Ba) = 0.015, /(^=^^Ba) = 0.107, fi^^'^Ba) = 0.080, /(^=^^Ba) = 0.758; 
and /(^^^Yb) = 0.180, /(^^^yb) = 0.219, /(^^=^Yb) = 0.191, /(^^^Yb) = 0.226, fi^'^^Yh) = 0.185 
We justify the s-process mixture choice below by showing that the total abundance distri- 
bution, involving 17 elements, follows an s-pro cess-dominant pattern. 

Pb I: The 4057.8 A line was easily detected in the syzygy spectrum (panel (d) of Fig- 
ure |8]), with EWo^tot ~ 5 mA. This line is very weak, but we confirmed its existence and 
approximate strength through co-addition of the six velocity-split spectra. The Pb abun- 
dance was derived from synthetic spectra in which the isotopic and hyperfine substructure 
were taken into account following Aoki et al. (2002), and using the Pb I line data of their 
Table 4. Variations in assumed Pb isotopic fractions produced no change in the derived 
elemental abundance. As with other features in the blue spectral region, CH contamination 
exists, but our syntheses suggested that it is only a small fraction of the Pb I strength. 

8. IMPLICATIONS OF THE CS 22964-161 ABUNDANCE PATTERN 

In the top panel of Figure [9] we plot the abundances relative to solar values for the 
10 n-capture elements detected in the CS 22964-161 syzygy spectrum, together with the 
n-capture abundances in the very s-process-rich BMP star CS 29497-030. These two stars 
have nearly the same overall [Fe/H] metallicity, as indicated by the horizontal lines in the 
panel. The n-capture overabundance pattern of CS 22964-161 clearly identifies it as another 
member of C- and s-process-rich "lead star" family. 

However, the relative n-capture abundance enhancements [X/Fe] of CS 29497-030 are 
about 1 dex higher than they are in CS 22964-161. This is emphasized by taking the dif- 
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ference between the two abundance sets, as is displayed in the lower panel of Figure [91 
These differences illustrate the trend toward weaker s-process enhancements of the heaviest 
elements in CS 22964-161 compared to CS 29497-030. 

Many studies have argued that these overall abundance anomalies, now seen in many low 
metallicity stars, must have originated from mass transfer from a (former) AGB companion 
to the stars observed today. Whole classes of stars with large C and s-process abundances 
are dominated by single-line spectroscopic binaries, including the high-metallicity "Ba II 
stars" (McClure, Fletcher, & Nemec 1980; McWilham 1988; McClure & Woodsworth 1990), 
the low-met aUicity red giant "CH stars" (McClure 1984), the "subgiant CH stars" (McClure 
1997), and the BMP stars (Preston & Sneden 2000; note that not all BMP stars share these 
abundance characteristics). The case is especially strong for the subgiant CH and BMP 
stars, for these are much too unevolved to have synthesized C and s-process elements in 
their interiors and dredged these fusion products to the surfaces. Our analysis has confirmed 
that primary and secondary CS 22964-161 stars are on or near the main sequence. Therefore 
we suggest that a third, higher-mass star is now or once was a member of the CS 22964-161 
system. During the AGB evolutionary phase of the third star, it transferred portions of its 
C- and s-pro cess-rich envelope to the stars we now observe. 

The very large Pb abundance of CS 22964-161 strengthens the AGB nucleosynthesis ar- 
gument. Gallino et al. (1998) and Travaglio et al. (2001) predicted substantial Pb production 
in s-process fusion zones of metal-poor AGB stars. In metal-poor stars, the neutron-to-seed 
ratio is quite high, permitting the neutron-capture process to run through to Pb and Bi, 
the heaviest stable elements along the s-process path. Prior to this theoretical prediction, it 
was assumed that a strong component of the s-process was required for the manufacture of 
half of the ^osp]^ ^-j^g solar system (Clayton & Rassbach 1967). That the patterns of the 
abundances of the neutron-capture elements in CS 22964-161 and CS 29497-030 resemble 
each other so well is a reflection of how easily low-metallicity AGB stars can synthesize the 
heavy elements. 

8.1. Inferred Nucleosynthetic and Dilution History 

We explored the origin of the neutron-capture enhancements in CS 22964-161 by com- 
paring the derived photospheric abundances with predicted stellar yields from the s-process. 
Employing FRANEC stellar evolutionary computations (see Straniero et al. 2003; Straniero, 
Gallino, & Cristallo 2006), we performed nucleosynthetic calculations following those of 
Gallino et al. (1998; 2006a) and Bisterzo et al. (2006). We then sought out good matches 
between the observed and predicted abundance pattern distributions with low-mass AGB 
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progenitors of comparable metallicities and a range of initial masses. 

We employed AGB models adopting different ^'^C-pocket efficiencies and initial masses 
to explore the nucleosynthetic history of the observed chemical compositions (Busso, Gallino 
& Wasserburg 1999; Straniero et al. 2003). Permitting some mixing of material between the 
s-process-rich contributions of the AGB donor, and the H-rich material of the atmosphere 
in which the contributions were deposited, we identified those AGB model progenitors that 
would yield acceptable fits for the predicted yields of all s-process elements beyond Sr. 

Ascertaining the best matches between the observed and predicted yields was performed 
in the following way. For a given initial AGB mass, we inspected the abundances of [Zr/Fe], 
[La/Fe], and [Pb/Fe] for various ^^C-pocket efficiencies adopted in the calculations. The 
difference between our predicted [La/Fe] and that of the observed abundance ratio gave us a 
first guess to the dilution of material for a given ^'^C-pocket efficiency. The dilution factor (dil) 
is defined as the logarithmic ratio of the mass of the envelope of the observed star polluted 
with AGB stellar winds and the AGB total mass transferred: dil = \og{Ml^^ / M^'^^^^). The 
^^C-pocket efficiency is defined in terms of the ^^C pocket mass involved in an AGB pulse 
that was adopted by Gallino et al. (1998). Their §2.2 states, "The mass of the ^^C pocket is 
5.0x10^^ Mq [here called ST], about 1/20 of the typical mass involved in a thermal pulse. 
It contains 2.8x10"^ Mq of ^^C." Thus our shorthand notation for ^'^C-pocket efficiency will 
be ST/N, where "N" is the reduction factor employed in generating a particular set of AGB 
abundance predictions. Subtracting this dilution amount from the predicted [Pb/Fe], we 
then compared the result to the observed [Pb/Fe]. The large range of ^^C-pocket efficiencies 
was then narrowed down, keeping only those results that fit the abundances of both [La/Fe] 
and [Pb/Fe] within 0.2 dex. Similar iterations were performed employing the abundances of 
[Zr/Fe]. 

We repeated this process for a range of initial AGB mass choices: 1.3, 1.5, 2, 3, and 
5M0. For AGB models of initial mass 3 and SM© no good match was found because the 
light s-elements (Is = Sr, Y, Zr) were predicted to have too high abundances with respect to 
the heavy s-elements (hs = Ba, La, Ce, Nd). For AGB models of M < 2Mq a satisfactory 
solution was found for Is, hs, and Pb, provided a proper ^'^C-pocket efficiency and dilution 
factor were chosen. In Figures [TU] and [TT] we show the matches between predicted and 
observed abundances for the two lowest-mass models, 1.3 and l.SM©, respectively. 

The abundances of the light elements Na and Mg further narrowed the range of allowable 
AGB progenitor models. Two independent channels are responsible for creating ^^Na: the 
^^Ne(p,7)^^Na reaction during H-shell burning, and the neutron capture on ^^Ne via the chain 
22Ne(ra,7)23Ne(/3-i/)23Na in the convective He flash (see Mowlavi 1999; Gallino et al. 2006). 
A large abundance of primary ^^Ne derives from the primary ^^C mixed with the envelope by 
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previous third dredge up episodes, then converted to primary ^"^N by HCNO burning in the H- 
burning shell and followed by double a capture via the chain ^'^N(a,7)^^F(/5+z/)^^0(a,7)^^Ne 
during the early development of the convective thermal pulse. The next third dredge up 
episode mixes part of this primary ^^Ne with the envelope. Finally, while the H-burning 
shell advances in mass, all the ^^Ne present in the H shell is converted to ^^Na by proton 
capture, accumulating in the upper region of He intershell. Note that in intermediate mass 
AGB stars suffering the so-called hot bottom burning (HBB) in the deeper layers of their con- 
vective envelope, efficient production of ^^Na results from proton capture on ^^Ne (Karakas 
& Lattanzio 2003). Furthermore, the marginal activation in the convective He flash of the 
reactions ^^Ne(Q;,n)^^Mg and ^^Ne(a,7)^^Mg would lead to enhanced Mg. At the same time, 
most neutrons released by ^^Ne(a,?T,)^^Mg are captured by the very abundant primary ^^Ne, 
thus producing ^^Na through the second channel indicated above. 

Comparison of the abundance data in Figures [TO] and [TT] shows that the observed [Na/Fe] 
argues for the exclusion of the l.SM© AGB model, which over-predicts this abundance ratio 
by ~0.5 dex. The same statement applies to the IMq model, not illustrated here. The 
larger Na abundances produced in these models, relative to the 1.3Mq model, is related to 
the larger number of thermal pulses (followed by third dredge up) that these stars experience. 
These arguments leave the I.BMq AGB model as the only one capable of providing a global 
good match to all observed elements. 

In order to match these abundances, the predicted AGB yields require an inferred dilu- 
tion by H-rich material. We do not know the distance between the AGB primary donor and 
the low mass companion (actually, the binary system) that is now observed. And, neither do 
we know precisely the mass loss rate of the AGB. Dilution by about 1 dex of s-process-rich 
AGB material with the original composition of the observed star, as we deduce on the basis 
of the above nucleosynthesis analysis, only fixes the ratio of the accreted matter with the 
outer envelope of the observed star. A more detailed discussion is taken up in §8.21 

A final remark concerns C, for which Figures [10] and [11] indicate an ~0.8 dex over- 
production compared to the observed [C/Fe]. This mismatch can in principle be substantially 
reduced by the occurrence of the so-called cool bottom process (CBB: see NoUett, Busso, 
& Wasserburg 2003; Dommguez et al. 2004a,b). This reduction would imply a substantial 
increase of predicted N abundance, from [N/Fe] ~ +0.5 shown in Figures [10] and [11] to 
perhaps +1.5. We attempted without success to detect the strongest CN absorption at the 
3883 A bandhead. Trial spectral syntheses of this wavelength region suggested that even 
[N/Fe] ~ +1.5 would not produce detectable CN. 
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8.2. Mixing and Dilution in the Stellar Outer Layers 

Very recently Stancliffe et al. (2007), extending the earlier discussion of Theuns, Boffin, 
& Jorissen (1996), report that extensive mixing due to thermohaline instability in a metal- 
poor star on the main sequence may severely dilute material that has been accreted from a 
companion. They computed the mixing time scale with the assumption that thermohaline 
convection behaved as a simple diffusion process, and they concluded that the new matter 
should rapidly be mixed down over about 90% of the stellar mass. 

Treating thermohaline convection simply as diffusion does not take into account the 
special nature of this process, which has extensively been studied in oceanography (e.g. 
Veronis 1965; Kato 1966; Turner 1973; Turner & Veronis 2000; Gargett & Ruddick 2003), 
and has also been considered in stars (Gough & Toomre 1982, Vauclair 2004). Thermohaline 
convection occurs in the ocean when warm salty water comes on top of cool fresh water. In 
this case the stabilizing thermal gradient acts against the destabilizing salt gradient. If the 
stabilizing effect compensates the destabilizing one, the medium should be stable, but it 
remains unstable due to double- diffusion: when a blob begins to fall down, heat diffuses 
out of it more rapidly than salt; then the blob goes on falling as the two effects no longer 
compensate. This creates the well known "salt fingers", thus is a very different physical 
environment than ordinary convection. 

A similar situation occurs in stars when high-/i matter comes upon a lower-/i region, 
which is the case for hydrogen-poor accretion. This has been studied in detail for the case of 
planetary accretion on solar- type stars (Vauclair 2004). As shown in laboratory experiments 
(Gargett & Ruddick 2003) the fingers develop in a special layer with a depth related to 
the velocity of the blobs and to the dissipation induced by hydrodynamical instabilities at 
their edges. This is a complicated process which may also be perturbed by other competing 
hydrodynamical effects. Depending on the situation, it is possible that the "finger" regime 
stops before complete mixing of the high-/x material into the low-/i one. In this case, we 
expect that the final amount of matter that remains in the thin sub-photospheric convective 
zone depends only on the final /i-gradient, whatever the original amount accreted: more 
accreted matter leads to more mixing so that the final result is the same. 

Another very important point has to be mentioned. At the epoch when accretion occurs 
on the dwarfs, these stars are already about 3 to 4 Gyr old. Gravitational settling of helium 
and heavy elements has already occurred and created a stabilizing /^-gradient below the 
convective zone (Vauclair 1999; Richard, Michaud, & Richer 2002). For example, as we 
illustrate in Figure [T2| in a 0.78Mq star with [Fe/H] = —2.3, after 3.75 Gyr the /i-value is as 
low as 0.58 inside the convective zone while it goes up to nearly 0.60 in deeper layers. Thus 
the stabilizing A/i//i is already of order 0.02. If the star accretes hydrogen-poor matter. 
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thermohaline convection may begin but it remains confined by the /x-barrier. The star can 
accrete as much high-/i matter as possible until this stabilizing /^-gradient is flattened. 

In §8.11 we saw, from the observed abundances and the AGB nucleosynthetic compu- 
tations, that the dilution factor is of the order of 1 dex for the AGB model of initial mass 
~1.3M0. Table [7] lists the predicted mass fraction of hydrogen (X), helium (Y), and heav- 
ier elements {Z) in the envelope of AGB progenitor stars of 1.3 and l.SM© for two cases. 
The first case describes the final abundance mix produced by the AGB at the end of its 
nucleosynthetic lifetime (the last Third Dredge- Up). The second case describes the mass 
average of the winds from the AGB over its s-processing lifetime. Also noted is the mean 
molecular weight of the material (/x), assuming fully ionized conditions. From Table [3, the 
/i-value in the wind lies in the range 0.61 to 0.78. From these computations, we find that, 
after dilution, the /i-value inside the convective zone should increase by A/i//i ~ 0.1% to 
4%. Clearly most of the inferred values for the dilution factor are compatible with stellar 
physics. For the lowest values it is possible that the accreted matter simply compensates the 
effect of gravitational diffusion. For larger values, thermohaline convection can have larger 
effects but, following Vauclair (2004), if we accept an inverse /x-gradient of order A/i//i = 
0.02 below the convective zone, such an accreted amount would still be possible. 



8.3. The Extraordinary Abundance of Lithium 

Near main-sequence-turnoff CEMP stars display a variety of Li abundances. The ma- 
jority have undetectable 6708 A Li I lines, implying log e(Li) < +1. A few have abundances 
similar to the Spite & Spite (1982) "plateau": log e(Li) = +2.10 ± 0.09 (Bonifacio et al. 
2007). These include our CEMP binary CS 22964-161, and at least CS 22898-027 (Thorburn 
& Beers 1992) and LP 706-7 (Norris, Ryan, & Beers 1997). A few more stars have Li abun- 
dances somewhat less than the Spite plateau value, e.g., log e(Li) ~ +1.7 in CS 31080-095 
and CS 29528-041 (Sivarani et al. 2006). This CEMP Li abundance variety stands in sharp 
contrast to the nearly constant value exhibited by almost all normal metal-poor stars near 
the main-sequence turnoff region. 

Although CS 22964-161 has a Spite plateau Li abundance, it is unlikely this is the same 
Li with which the system was born. Considering just the primary star, suppose that the C 
and s-process-rich material that it received was Li-free. If this material simply blanketed the 
surface of the primary, mixing only with the original Li-rich atmosphere and outer envelope 
skin, then the resulting observed Li would be at least less than the plateau value. Alter- 
natively, if the transferred material induced more significant mixing of the stellar envelope, 
the surface Li would be severely diluted below the limit of detectability in the primary's 
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atmosphere. In order to produce the large observed Li abundance in the primary, freshly 
minted Li must have been transferred from the AGB donor star. 

Li can be produced during the AGB phase both in massive AGB stars and in low-mass 
AGB stars via ^He(a,7)^Be(e~, i^e) ''Li, but hot protons rapidly destroy it via ''Li(p,a)'^He. 
Cameron & Fowler 1971) proposed that in some circumstances the freshly minted ^Be can 
however be transported to cooler interior regions before decaying to ^Li. In massive AGB 
stars ( 4:Mq ^ M < 7Mq) nucleosynthesis at the base of the convective envelope ("hot 
bottom burning" ; Scalo, Despain, & Ulrich 1975) can produce abundances via the Cameron- 
Fowler mechanism as high as log e(^Li) ~ 4.5 (Sackmann & Boothroyd 1992). However, the 
lifetimes of these massive AGB stars are far too short to be considered as a likely source of 
Li in CS 22964-161. 

In models of low-mass AGB stars {IMq < M < 3Mq) of [Fe/H] = -2.7, Iwamoto et al. 
(2004) find that a H-flash episode can occur subsequent the first fully developed He shell 
flash. The convection produced by the He shell flash in this metal-poor star model can reach 
the bottom of the over-lying H-layer and bring protons down into the He intershell region, a 
region hot enough to induce a H flash. The high-temperature conditions of the H-fiash can 
then produce ''Li by the Cameron-Fowler mechanism. Although the Iwamoto et al. 2.5M0 
model produced a weak H flash, the less massive models produced more energetic ones. In 
the lowest mass models, surface abundances of log e(''Li) > 3.2 were achieved. Note that 
there is a maximum metallicity close to [Fe/H] = —2.5 beyond which such an H-flash is not 
effective at producing ^Li. Another viable mechanism for producing a very high abundance 
of ^Li in low mass AGB stars of higher metallicities, as observed in the intrinsic C(N) star 
Draco 461 ([Fe/H] = —2.0 ± 0.2), was advanced by Dommguez (2004a,b) and is based on 
the operation of the CBP introduced in §8.1[ 

Although s-process calculations were not performed in the Iwamoto et al. study, we note 
that it would seem that the parameters of the best-fit model of §8. II would produce sufficient 
Li to fit the observations of CS 22964-161. 



9. IMPLICATIONS FOR THE PUTATIVE AGB DONOR STAR 

The binary nature of CS 22964-161 presents a unique opportunity to explore the nature 
of the AGB star presumably responsible for its abundance anomalies. The ingredients of 
the exploration are: (a) the analysis of our radial velocities of CS 22964-161 obtained over 
an interval of about 1100 days, (b) estimation of the orbit changes caused by accretion 
from an AGB companion, (c) stability of the hierarchical triple system in its initial and 
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final states (Donnison & Mikulskis 1992, hereafter DM92; Szebeliely & Zare 1977; Tokovinin 
1997; Kiseleva, Eggleton, & Anolova 1994), and (d) application of an appropriate initial-final 
mass relation for stars of intermediate mass (Weidemann 2000, and references therein). 

9.1. Evidence for an AGB relic? 

Velocity residuals calculated for our orbit solution show no trend with time over the 
nearly 1100-day interval of observation. This is illustrated by the plot of orbital velocity 
residuals versus Julian Date in Figure [131 Thus, we have yet to find evidence for a third 
component in the system. The Julian Date interval of our CS 22964-161 data is modest 
compared to the longest orbital periods of the so-called Barium and CH giant stars (McClure 
& Woodsworth 1990) and main-sequence/subgiant CH stars (McClure 1997) summarized in 
Table [61 On the other hand the errors of our observations (cr = 0.5 km s~^) are small 
compared to typical velocity semi-amplitudes of these comparison C-rich stellar samples, so 
drift of the center-of-mass velocity by as much as 1 km s~^ due to the presence of a third star 
should be detectable. Accordingly, Figure [131 encourages us to contemplate the possibility 
that CS 22964-161 is not accompanied by a white dwarf relic of AGB evolution. Continuing 
observation will be required to test this notion. 

9.2. Change in Orbit Dimensions Induced by Accretion 

The consequences of mass accretion by a binary were first explored by Huang (1956), who 
calculated the change in semi-major axis of a binary embedded in a stationary interstellar 
cloud. Using McCrea's (1953) "retarding force" that followed from Bondi-Hoyle-Lyttleton 
accretion theory (Bondi & Hoyle 1944), Huang found that mass accretion shrinks orbits, and 
he proposed that this process could produce close binaries. The problem of binary accretion 
has been revisited in recent years by numerical simulations of star formation in molecular 
clouds. Bate & Bonnell (1997) calculate the changes in binary separation due to accretion 
by an infalling circumstellar cloud. The change is very sensitive to the assumed specific 
angular momentum (SAM) carried by the accreted material. In the case of low SAM they 
recover Huang's result. For large SAM orbits the binaries expand. More recently Soker 
(2004) has calculated the SAM of mass accreted by a binary system embedded in the wind 
of an AGB star, which is our case. For triple systems containing a distant AGB companion 
with approximately coplanar inner and outer orbital orbits (statistically most probable in 
the absence of a special orientation mechanism) he also recovers Huang's result. Finally, 
building on an earlier analysis by Davies & Pringle (1980), Livio et al. (1986) found that the 
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rate of accretion of angular momentum from winds is significantly less than that deposited 
in the simple Bondi-Hoyle accretion model. 

Our gloss on all these analyses is that the most likely consequence of mass transfer to 
CS 22964-161 from its AGB companion a long time ago was a decrease in separation, i.e., 
the initial separation was larger than its present value of ~0.9 AU. Conservatively, we use 
0.9 AU as a lower limit in the estimates that follow below. 

Empirical stability criteria for hierarchical triple systems are expressed as the ratio of 
the semi-major axis of the outer orbit to the semi-major axis of inner orbit (Heintz 1978; 
Szebehely & Zare 1977), and alternatively as the ratio of periods of the outer and inner 
orbits (Kiseleva et al. 1994; Tokovinin 1997). If these ratios exceed critical values, a triple 
system can be long-lived. Note, however, the cautionary remarks of Szebehely & Zare about 
the role of component masses in evaluation of stability for particular cases, and the relevance 
of these masses to the conclusions of DM92 about the fate of our AGB relic. 

Let ai2 be the semi-major axis of the relative orbit of CS 22964-161 and 03 be the relative 
semi-major axis that joins the center of mass of CS 22964-161 to the S^'^ (AGB-to-be) star. 
If we accept the Heintz (1978) and Szebehely & Zare (1977) stability criterion a-^jayi > 8, 
then as > 7.3 AU and the orbital period of the outer binary had to be P > 4000 d. For 
Tokovinin's (1997) more permissive criterion, P3/P12 > 10 the minimum period is P > 2500 d 
and 03/012 > 4.6. Under these conditions none of our candidate metal-poor AGB stars (see 
the discussion in ^KB with M < 2Mq fill the Roche lobes (Marigo & Girardi 2007) of 
these minimum orbits, so wind accretion must have been the mechanism of mass transfer. 
According to Theuns et al. (1996) a 3Mq AGB star in a 3-AU circular orbit will transfer no 
more than ~1% of its wind mass to a 1.5 solar mass companion by wind accretion. Applying 
their result to our case, M(AGB) < 1.5Mq and AGB relic mass ~O.5M0, the maximum 
mass transfer is about (1.5— 0.5) xO. 015 ~ O.O15M0. 

The remainder of the AGB wind mass is ejected from the system. We calculated the 
change in period and semi-major axis of the AGB orbit following Hilditch (2001) for the case 
of a spherically symmetric wind. The orbit expands and the period lengthens according to 

ln(a//aj) = ln(l + Ui) — ln(l + uj) 

and 

ln(P//Pi) = 2 X ln(l + Mi) - 2 X ln(l + Uf) 

in which initial and final values of u are Ui = M^i/lMi + M2) and Uf = M^f/lMi + M2) . 
We have ignored small changes in Mi and M2. For M3 = 1.3 Mq and the Tokovinin (1997) 
stability condition (03/012 > 4.6), the minimum orbit expands from 4.2 AU to 5.9 AU and 
the minimum period lengthens from 1900 days to 3700 days. For the more restrictive Heintz 
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{a^/ai2 > 8) stability condition the minimum orbit expands from 7.3 AU to 10.2 AU and the 
minimum period lengthens from 4300 days to 8400 days. The K- values for the center-of-mass 
of the CS 22964-161 binary for the two cases are 4.3 km s~^ and 3.3 km s~^, respectively, 
both readily detectable by conventional echelle spectroscopy. 

9.3. Dynamical Stability of CS 22964-161 as a Hierarchical Triple System 

The DM92 study investigated the stability of triple systems that consist of a close inner 
binary with masses (Mi,M2) and semi-major axis al attended by a remote companion of 
mass M3 at distance q2 from the center of mass of the inner binary. They use the parameter 
notation of Harrington's (1975) pioneering numerical simulations. DM92 investigate stability 
in three mass regimes, by integration of the equations of motion for various sets of masses, 
initial positions and velocities: 

case (i): M3 < minimum ( Mi, M2) 

case (ii): Mi > M3 > M2 

case (Hi): M3 > maximum ( Mi, M2) 
They pursued the integrations for at least 1000 orbits of the inner binary or until disruption, 
which ever happened first. Among the several results in §5 of their paper one is of particular 
interest for us: in their words "When M3 was the least massive body in the system, the 
system invariably became unstable through the tendency of M3 to escape from the system 
altogether." Accepting the McClure (1984) paradigm, we suppose that CS 22964-161 began 
its existence (in) hierarchical triple system in which a relatively low-mass AGB 

progenitor ( §8.11) was orbited by a close binary with component masses only slightly smaller 
than their present values, Misin^i = 0.77Mq and M2sin^i = 0.68 Mq. Evidently, the distance 
between M3 and the close binary was large enough to insure stability during the lifetime 
of M3. Curiously, although Kiseleva et al. (1994) made numerical simulations apparently 
similar to those of DM92, they make no mention of the ejection of low-mass outer components 
that attracted our attention to DM92. 

We confine our attention to initial values of AGB mass <2.OM0, adopted in accordance 
with the nucleosynthetic calculations in §8.1[ In the late stages of AGB evolution the system 
evolved from case (Hi) to case (i) via mass loss due to the AGB superwind. For all acceptable 
models in §8. li the mass of the WD relic for an AGB initial mass model of I.3M0 is less than 
O.6OM0 (Weidemann 2000). Therefore, according to DM92 M3 will be ejected following this 
ancient AGB evolution. Our extant radial velocity data suggest that the ejection already 
occurred. From the effect of the AGB relic on the 7 velocity of CS 22964-161 calculated 
above we believe that this conclusion can be tested by future observations. 
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10. CONCLUSIONS 

We have obtained high resolution, high S/N spectra over a three- year period for CS 22964-161, 
which was known to be metal-poor and C-rich from previous low resolution spectroscopic 
studies. We discovered CS 22964-161 to be a double-line spectroscopic binary, and have de- 
rived the binary orbital parameters (P = 252 d, e =0.66, Mp sin^ i = 0.77 Mq, Mg sin^ i = 
0.68 Mq). Both binary members lie near the metal-poor main-sequence turnoff {T^s^p = 
6050 K, log gp = 3.7, and Tcs,s =5850 K, log gs = 4.1). We derived similar overall metal- 
licities for primary and secondary, [Fe/H] = —2.4 and similar abundance ratios [X/Fe]. In 
particular, both stars are similarly enriched in C and the n-capture elements, with a clear 
s-process nucleosynthesis signature. The primary has a large Li content; the secondary's Li I 
feature is undetectably weak, and does not usefully constrain its Li abundance. 

The observed Li, C and s-process abundances of the CS 22964-161 system must have 
been produced by an AGB star whose relic, depending on its mass, may have been ejected 
from the post-AGB hierarchical triple. It seems AGB enrichment was produced by minor 
mass accretion in the AGB super-wind rather than Roche-lobe overflow. As discussed in 

and contrary to the thermohaline diffusion suggestion of Stancliffe et al. (2007), such 
hydrogen-poor accreted matter probably remained in the outer stellar layers, owing to the 
stabilizing /i-gradient induced by helium gravitational settling. Strong thermohaline diffu- 
sion is difficult to reconcile with the observed Li: about 9/10 of the stellar mass would be 
efficiently mixed on a time scale of 1 Gyr, which would imply the complete destruction of ^Li 
because of the very high temperature reached in the inner zones. It seems more likely that 
the moderate thermohaline mixing discussed ^ related to the effect of gravitational settling 
in the first 3 to 4 Gyr before mass accretion by the AGB donor, would save the Li from de- 
struction. In the model of a low mass AGB of low metallicity calculated by Domfnguez et al. 
(2004a,b) with the inclusion of "cool bottom processing" during the AGB phase, log e(Li) > 
3 is achieved in the envelope. A dilution factor of about 1 dex would bring the Li abundance 
into accord with the observed value. 
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Fig. 1. — Top panel: radial velocity observations of CS 22964-161 primary and secondary 
stars, and best-fit orbital solutions from these velocities. Circles represent data taken with 
MIKE on the Magellan Clay telescope and squares represent data taken with the echelle 
spectrograph on the du Pont telescope. Filled symbols represent data for the primary and 
open symbols are for the secondary. The lines are fits to the Magellan data only. Bottom 
panel: differences of the observed velocities and the orbital solutions. 
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Fig. 2. — Observed composite absorption spectra of primary and secondary stars. The 
bottom spectrum is the spectrum at syzygy (AVr = km s"-*^); this is the mean of three 
individual observations. The middle spectrum is the mean of three observations obtained at 
phases with AVr = 60 km s~^. The top spectrum represents the original syzygy spectrum 
after application of a velocity shift of 60 km and dilution by addition of a continuum flux 
that yields secondary line depths approximately matching those seen in the middle spectrum. 
Dotted lines point to secondary positions of a few lines that are very strong in the primary 
spectrum. The flux scale of the syzygy spectrum is correct, and that other two spectra have 
been vertically shifted in the figure for display purposes. The rest wavelength scale is that 
of the primary star. 
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Fig. 3. — Comparison of observed equivalent widths of the secondary star inferred from the 
syzygy spectrum and those measured directly from the six spectra with large velocity splits 
between primary and secondary stars. Measurements for each spectral feature are connected 
by vertical hues and labeled with the feature wavelength. Medians of the EW% are displayed 
as X symbols. The dashed slanting line indicates equality of EW values. 



-34- 



6, 
5, 
5, 
4, 

^ 4, 
o " 

5, 
5, 
4, 
4, 



- 



5 



5 




5 



5 





|_ o o o 

00 



o Fe I 
X Fe II 



O ■ ■ ■ ■ n- ■ ■ ■ 



primary 



± 



o 

o ^ 
o 



o o 
o 



.o. 



OOqO 

ho-^" o 

X o 
o ^ o 



± 



o o 
o o 



cx> 



^ o o 

o o 
o O 



secondary 



± 



4000 5000 

Wavelength (A) 



6000 



Fig. 4. — Iron abundances plotted as a function of wavelength for primary (top panel) 
and secondary (bottom panel) stars. The symbols are defined in the figure legend. Dashed 
horizontal lines are drawn to indicate the mean Fe I abundances of primary and secondary 
stars, log e(Fe) = 5.13 and 5.11, respectively. 
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Fig. 5. — Evolutionary tracks (VandenBerg et al. 2006) and CEMP data plotted in the 
log Tefj — log g plane near the metal-poor main-sequence turnoff. The chosen tracks corre- 
spond to three different inclinations of the binary; heavy and light lines are for the primary 
and secondary, respectively. A full track covers an age range of 4 to 14 Gyr. The figure leg- 
end identifies the masses of the tracks for three different CS 22964-161 orbital sin i choices, 
using the notation (sm i: Mp, Mg). Filled circles represent the CS 22964-161 primary and 
secondary (log Tgfr, log g) values derived from the spectroscopic analysis ( §5.ip . The plotted 
error bars for the primary are also from the spectroscopy. The square shows the secondary 
parameters from the dimensional analysis described in §5.21 The line about this point corre- 
sponds to ATofj = ±400 K. The triangles represent CEMP stars from the literature; see §5.21 
for details. The triangles for three pairs of stars with identical temperatures and gravities 
have been shifted apart from each other for display purposes. 
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Fig. 6. — Observed and synthetic spectra of the Sr II 4077.71 A hne in observation 13817 
of CS 22964-161 (AVr ~ 58.09 km s"^). The top spectrum shows the full relative flux scale, 
and the bottom spectrum covers just the relative flux region above 0.75. The wavelength 
scale is at rest velocity for the primary star. The observed spectrum is depicted with open 
circles. The blue line represents a synthetic spectrum computed without any contribution 
from Sr. The black line shows the best overall fit to primary and secondary Sr II features: 
log e(Sr) = 0.98 and 0.78, respectively. The red and green lines indicate syntheses for each 
star that are 0.5 dex smaller and larger, respectively, than the best fits. 
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Fig. 7. — Observed and synthetic spectra of the CH G-bandhead in CS 22964-161. The 
observation is 13817, for which AVr = 58.09 km s~^. Open circles represent the observed 
spectrum. In the top panel only the contribution of the primary star to this blend is con- 
sidered. The abundances of the synthetic spectra, in order of increasing CH strength, are 
log e(C)p = — oo, 7.24, 7.64, and 8.04 (the blue, red, black, and green lines, respectively). 
In the bottom panel, the abundance of the primary is fixed at log e(C)p = 7.64, and that of 
the secondary is, again in order of increasing CH strength, log e(C)s = — oo, 6.44, 7.44, and 
8.44. An arrow indicates the wavelength offset between primary and secondary stars of this 
observation. 
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Fig. 8. — Selected spectral features of n-capture species in the CS 22964-161 syzygy spec- 
trum. 
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Fig. 9. — In the top panel, n-capture abundances relative to solar values [X/H] from the 
CS 22964-161 combined spectrum are plotted along with those of the BMP very s-process- 
rich star CS 29497-030 (Ivans et al. 2005). Horizontal lines are drawn to indicate the general 
Fe-metallicity levels of the two stars. In the bottom panel the abundance differences between 
these two stars are shown. 
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Fig. 10. — Comparison of the [X/Fe] abundances in CS 22964-161 with predictions from 
s-process calculations of a I.SMq AGB star model. The solid red line corresponds to the best 
match between the observed and predicted abundance pattern. The dotted blue and dashed 
green lines indicate the difference that adopting models of iO.OSM© would make (effectively, 
changing the number of thermal pulses from 4 to 6, identified in the figure legend as n4, n5, 
n6). For the specific choice of the ^^C-pocket (ST/12), and for the definition of the dilution 
factor (dil), see the text. 
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Fig. 11. — Comparison of the [X/Fe] abundances in CS 22964-161 with predictions from 
s-process calculations of a 1.5Mq AGB star model. The solid red line corresponds to the 
best match between the observed and predicted abundance pattern. The dotted magenta, 
long-dashed green, and short-dashed blue lines show the difference made to the pattern 
by adopting different amounts of dilution. For discussion of the the various choices of the 
^^C-pocket efficiency (ST/N), and the dilution factor (dil), see the text. 
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Fig. 12. — Mean molecular weight (/x) and temperature (T) profiles, as functions of the 
fractional mass (mass above the considered layer), in a O.TSM© star with [Fe/H] = -2.3, at 
an age of 3.75 Gyr. Here /i is the real molecular weight, including partial ionization. The 
fractional mass of the convective zone is 0.004. The helium depletion due to gravitational 
settling, of about 20% in this case, already leads to an important stabilizing yU-gradient 
below the convective zone. This model is one of those computed by Richard ct al. (2002), 
including pure atomic diffusion. The full model was made available to us by Richard (private 
communication) . 
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Fig. 13. — Velocity residuals to the fitted orbit are plotted against the date of observation 
for the Magellan data. Filled circles represent the primary and open circles represent the 
secondary. These residuals show no trend with time. 
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Table 1. Velocity Measurements for CS22964-161 



HJD 




v. 


Phase 


(- 2440000) 


km s~^ 


km s~^ 




du Pont Observations 


13592.52532 


24.79 


41.04 


0.0716 


13592.56468 


24.63 


42.35 


0.0718 


13593.52234 


25.34 


38.32 


0.0756 


13593.56091 


25.21 


40.26 


0.0757 


13867.92873 


36.59 


27.00 


0.1623 


13870.91114 


36.44 


25.90 


0.1741 


13871.85817 


37.22 


26.27 


0.1778 


Magellan Observations 


13230.66465 


38.43 


25.01 


0.6386 


13580.79663 


5.09 


65.21 


0.0252 


13585.74591 


14.80 


52.79 


0.0448 


13587.62635 


18.56 


49.77 


0.0522 


13602.62780 


32.29 




0.1116 


13603.50474 


33.00 




0.1151 


13606.50536 


33.83 




0.1270 


13631.67197 


39.08 


26.29 


0.2266 


13634.54878 


39.06 


25.43 


0.2380 


13817.90625 


6.10 


64.19 


0.9642 


13832.88449 


4.11 


66.51 


0.0235 


13834.88836 


8.21 


60.76 


0.0314 


13889.88563 


40.09 


25.83 


0.2492 


13890.73000 


40.19 


26.24 


0.2526 


13891.72338 


40.03 


25.80 


0.2565 


13891.79617 


39.87 


25.63 


0.2568 


13891.91199 


39.79 


25.33 


0.2572 


13892.72313 


40.47 


25.77 


0.2604 


13892.91348 


40.03 


25.21 


0.2612 


13893.76837 


40.02 


25.34 


0.2646 


13897.80466 


40.22 


24.87 


0.2806 
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Table 1 — Continued 



HJD Vp V, Phase 

(- 2440000) km s"^ km s"^ 



13898. 


,72511 


40, 


.00 


24, 


.35 


0, 


.2842 


13917. 


,74748 


40, 


.82 


23, 


.51 


0, 


.3595 


13919. 


.72723 


40, 


.20 


23 


.00 


0, 


.3674 


13935. 


.76042 


41, 


.08 


23, 


.66 


0, 


.4309 


13938. 


.71168 


41, 


.14 


23, 


.29 


0, 


.4426 


13938. 


.73591 


41, 


.10 


23 


.53 


0, 


.4427 


13951. 


,57253 


40, 


.34 


23 


.09 


0, 


.4935 


13989. 


,60930 


38, 


.73 


25, 


.04 


0, 


.6441 


13990. 


,59320 


38, 


.69 


25, 


.36 


0, 


.6480 


14197. 


,85590 


40, 


.80 


23 


.87 


0, 


.4689 


14258. 


,82674 


37, 


.60 


28, 


.04 


0, 


.7088 


14259. 


.78553 


37, 


.50 


27, 


.81 


0, 


.7126 


14316. 


.78097 


15, 


.19 


52, 


.81 


0, 


.9383 


14317. 


.65762 


13, 


.63 


53, 


.75 


0, 


.9418 


14328. 


,67079 


-4, 


.33 


74, 


.49 





.9854 


14329. 


,70221 


-5, 


.60 


75, 


.71 


0, 


.9895 


14333. 


,56647 


-5, 


.03 


75, 


.31 


0, 


.0048 


14334. 


,49082 


-3 


.13 


74, 


.78 


0, 


.0085 


14335. 


,48607 


-1, 


.45 


72, 


.22 


0, 


.0124 


14336. 


.53932 


0, 


.46 


69, 


.88 


0, 


.0166 
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Table 2. Orbital Parameters for CS22964-161 



Parameter 


Value 


P (days) 


252.481±0.043 


To (HJD-244 0000) 


13827.383±0.103 


7 {km s~^) 


32.85± 0.05 


Kp {km s^^) 


23.65±0.10 


Kg {km s~^) 


26.92±0.13 


e 


0.6564±0.0024 


uj (deg) 


188.54±0.39 


CTp {km s~^) 


0.37 


(Tg {km s~^) 


0.56 


Derived quantities: 




A sin i {Rq) 


190.32±0.65 


Mp sin^ i {Mq) 


0.773±0.009 


Ms sin^ i (Mq) 


0.680±0.007 



-47- 



Table 3. Line List for CS 22964-161 



Species A x log gf EWp^o EWp,t EWs,o EWs,t 

A eV mA mA mA mA 



Li I 


6707.80 


0, 


.00 


+0 


.18 


24 


29 






Mg I 


3829.36 


2, 


.71 


— 0, 


.21 


104 


122 


31 


222 


Mg I 


3832.31 


2, 


.71 


+0, 


.14 


128 


150 


30 


213 


Mg I 


4703.00 


4, 


.35 


— 0. 


.38 


40 


47 


8 


55 


Mg I 


5172.70 


2, 


.71 


— 0. 


.38 


127 


152 


25 


156 


Mg I 


5183.62 


2 


.72 


— 0. 


.16 


144 


172 


28 


174 


Mg I 


5528.42 


4 


.35 


— 


.34 


40 


49 


10 


63 


Al I 


3961.53 


0, 


.01 


— 0, 


.34 


62 


73 


11 


79 


Si I 


3905.53 


1, 


.91 


+ 


.09 










Ca I 


4226.74 


0. 


.00 


+0, 


.24 


134 


158 


27 


180 


Ca I 


5349.47 


2, 


.71 


— 0. 


.31 


11 


13 






Ca I 


5588.76 


2, 


.53 


+0, 


.36 


21 


28 


3 


27 


Ca I 


5857.46 


2, 


.93 


+0, 


.24 


17 


21 






Ca I 


6102.73 


1, 


.88 


— 0, 


.79 


12 


14 


4 


23 


Ca I 


6122.23 


1, 


.89 


— 0, 


.32 


26 


33 


2 


12 


Ca I 


6162.18 


1, 


,90 


— 0, 


.09 


33 


41 


4 


27 


Ca I 


6439.08 


2, 


.52 


+0, 


.39 


28 


36 


3 


22 


Sc II 


3645.31 





,02 


— 


.42 


43 


50 


8 


61 


Sc II 


4246.84 





.31 


+0, 


.24 


59 


69 


8 


55 


Ti II 


3477.19 


0, 


.12 


— 1, 


.06 


85 


99 






Ti II 


3510.85 


1, 


.89 


+0, 


.14 


55 


64 






Ti II 


3596.05 


0, 


.61 


— 1, 


,22 


43 


50 


12 


90 


Ti II 


3641.34 


1, 


.24 


— 0, 


.71 


40 


46 


12 


86 


Ti II 


3759.30 


0, 


.61 


+0, 


,20 


97 


113 


19 


137 


Ti II 


4394.06 


1, 


.22 


— 1, 


,77 


16 


18 






Ti II 


4443.80 


1, 


.08 


— 0, 


,70 


55 


66 


6 


45 


Ti II 


4444.56 


1, 


.12 


— 2, 


.21 


14 


16 






Ti II 


4468.49 


1, 


.13 


— 0, 


.60 


59 


71 


7 


52 


Ti II 


4501.27 


1, 


.12 


— 0, 


.76 


54 


65 


5 


40 


Ti II 


5336.79 


1, 


.58 


— 1, 


.63 


13 


16 






Cr I 


3578.69 


0, 


.00 


+0, 


,41 


59 


69 


13 


97 


Cr I 


3593.50 


0. 


.00 


+0, 


,31 


65 


76 


6 


46 


Cr I 


4254.33 


0. 


.00 


— 0, 


,11 


57 


67 


4 


29 


Cr I 


4274.80 


0, 


.00 


— 0, 


.23 


72 


84 






Cr I 


4289.72 


0, 


,00 


— 0, 


.36 


63 


74 


7 


44 


Cr I 


4646. 17 


1, 


.03 


— 0, 


.73 


5 


6 


2 


14 


Cr I 


5206.04 


0, 


.94 


+ 


,03 


29 


35 


6 


36 


Cr I 


5409.79 


1, 


.03 


— 0, 


,71 


15 


18 






Cr II 


4558.65 


4, 


.07 


— 0, 


,66 


8 


9 


2 


10 


Cr II 


4588.20 


4, 


.07 


— 0, 


,64 


6 


7 






Mn I 


4030.76 


0, 


.00 


— 0, 


.48 










Mn I 


4033.06 


0, 


.00 


— 0, 


.62 










Mn I 


4034.49 


0, 


.00 


-0 


.81 










Fe I 


3475.46 


0, 


.09 


-1 


.05 


82 


94 


34 


249 


Fe I 


3476.71 


0, 


.12 


-1, 


.51 


73 


84 


8 


61 


Fe I 


3490.59 


0, 


.05 


-1, 


,11 


76 


88 


16 


117 


Fe I 


3497.84 


0, 


.11 


-1, 


,55 


71 


82 


8 


58 


Fe I 


3521.27 


0, 


.91 


-0, 


.99 


56 


65 


15 


111 


Fe I 


3554.94 


2, 


.83 


+ 


,54 


47 


55 


6 


45 


Fe I 


3558.53 


0, 


,99 


-0, 


,63 


68 


79 


16 


113 


Fe I 


3565.40 


0, 


,96 


-0, 


,13 


89 


103 


21 


151 


Fe I 


3570.13 


0, 


.91 


+0. 


,18 


151 


175 


26 


190 


Fe I 


3581.21 


0. 


.86 


+0, 


.42 


114 


132 


36 


262 


Fe I 


3606.69 


2, 


.69 


+0, 


,32 


41 


47 


8 


59 


Fe I 


3608.87 


1, 


.01 


-0, 


.09 


80 


92 


20 


142 


Fe I 


3618.78 


0, 


.99 


+ 


,00 


88 


102 


30 


219 


Fe I 


3631.48 


0, 


.96 


+0, 


,00 


102 


118 


23 


168 


Fe I 


3647.85 


0, 


,91 


-0, 


,14 


82 


95 


22 


162 


Fe I 


3679.92 


0, 


.00 


-1, 


,58 


69 


80 


17 


125 


Fe I 


3687.47 


0, 


.86 


-0, 


,80 


87 


102 


11 


77 


Fe I 


3727.63 


0, 


.96 


-0, 


,61 


83 


96 


10 


75 


Fe I 


3745.57 


0, 


.09 


-0, 


.77 


100 


116 


25 


176 


Fe I 


3745.91 


0, 


.12 


-1, 


.34 


89 


104 


10 


72 


Fe I 


3758.24 


0, 


.96 


-0, 


.01 


100 


116 


21 


149 


Fe I 


3763.80 


0, 


.99 


-0, 


,22 


86 


99 


15 


106 
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Table 3 — Continued 



Species 


X 


X 


log gf 


EWp^a 






EWs,t 




A 


cV 




mA 


mA 


mA 


mA 



Fe I 


3787.89 


1 


01 


-0 


84 


69 


80 


10 


71 


Fe I 


3820.44 





86 


+0 


16 


114 


133 


22 


157 


Fe I 


3826.89 





92 


-0 


02 


99 


116 


22 


153 


Fe I 


3840.45 





99 


-0 


50 


74 


87 


20 


140 


Fe I 


3856.38 





05 


-1 


28 


84 


98 


11 


76 


Fe I 


3865.53 


1 


01 


-0 


95 


67 


79 


11 


76 


Fe I 


3899.72 





09 


-1 


51 


70 


81 


19 


135 


Fe I 


3917.18 





99 


-2 


15 


27 


31 


4 


27 


Fe I 


3922.92 





05 


-1 


63 


76 


87 


16 


110 


Fe I 


3949.96 


2 


18 


-1 


26 


18 


21 


6 


38 


Fe I 


4005.25 


1 


56 


-0 


58 


61 


71 


12 


85 


Fe I 


4071.75 


1 


61 


-0 


01 


79 


92 


17 


125 


Fe I 


4134.69 


2 


83 


-0 


65 


19 


22 


4 


29 


Fe I 


4143.88 


1 


66 


-0 


61 


68 


80 


8 


53 


Fe I 


4147.68 


1 


49 


-2 


07 


12 


14 


4 


29 


Fe I 


4175.64 


2 


85 


-0 


83 


18 


21 






Fe I 


4187.05 


2 


46 


-0 


61 


38 


46 


2 


16 


Fe I 


4187.81 


2 


42 


-0 


51 


42 


49 


6 


41 


Fe I 


4191.44 


2 


47 


-0 


67 


26 


30 


8 


54 


Fe I 


4199.10 


3 


05 


+0 


16 


52 


61 






Fe I 


4202.04 


1 


49 


-0 


69 


62 


72 


13 


85 


Fe I 


4216.19 





00 


-3 


36 


24 


28 


2 


13 


Fe I 


4222.22 


2 


46 


-0 


91 


24 


28 






Fe I 


4233.61 


2 


48 


-0 


68 


37 


43 


4 


24 


Fe I 


4250.13 


2 


47 


-0 


38 


47 


56 






Fe I 


4260.48 


2 


40 


+0 


08 


61 


72 


10 


69 


Fe I 


4271.16 


2 


45 


-0 


34 


58 


69 


20 


131 


Fe I 


4271.77 


1 


49 


-0 


17 


93 


109 


11 


77 


Fe I 


4282.40 


2 


18 


-0 


78 


34 


39 


8 


52 


Fe I 


4383.56 


1 


48 


+0 


21 


104 


122 


16 


100 


Fe I 


4404.76 


1 


66 


-0 


16 


76 


90 


14 


92 


Fe I 


4415.12 


1 


61 


-0 


62 


62 


73 


16 


108 


Fe I 


4427.32 





05 


-2 


92 


37 


44 


10 


65 


Fe I 


4430.61 


2 


22 


-1 


73 


10 


12 


2 


15 


Fe I 


4447.72 


2 


22 


-1 


34 


19 


22 


4 


27 


Fe I 


4461.66 





09 


-3 


19 


20 


24 


8 


52 


Fe I 


4466.55 


2 


83 


-0 


60 


19 


22 


9 


59 


Fe I 


4489.73 





12 


-3 


90 


6 


6 


3 


20 


Fe I 


4494.56 


2 


20 


-1 


14 


23 


27 


3 


20 


Fe I 


4628.61 


2 


18 


-0 


89 


35 


42 


8 


51 


Fe I 


4602.96 


1 


49 


-2 


21 


13 


15 


2 


16 


Fe I 


4871.32 


2 


87 


-0 


36 


26 


29 


5 


34 


Fe I 


4872.14 


2 


88 


-0 


67 


26 


30 






Fe I 


4890.76 


2 


88 


-0 


39 


24 


29 


10 


65 


Fe I 


4891.50 


2 


85 


-0 


11 


37 


44 


8 


49 


Fe I 


4903.32 


2 


88 


-0 


93 


12 


14 






Fe I 


4919.00 


2 


87 


-0 


34 


29 


34 


6 


35 


Fe I 


6006.12 


2 


83 


-0 


61 


23 


28 






Fe I 


6014.94 


3 


94 


-0 


27 


12 


16 


4 


22 


Fe I 


6049.83 


2 


28 


-1 


36 


18 


21 


8 


49 


Fe I 


5133.69 


4 


18 


+0 


20 


11 


13 


7 


45 


Fe I 


5191.47 


3 


04 


-0 




21 


25 






Fe I 


5192.35 


3 


00 


-0 


42 


30 


36 






Fe I 


5194.95 


1 


56 


-2 


02 


17 


25 






Fe I 


6202.34 


2 


18 


-1 


84 


11 


13 






Fe I 


6216.28 


1 


61 


-2 


08 


16 


19 






Fe I 


6232.96 


2 


94 


-0 


06 


39 


48 


3 


29 


Fe I 


6266.66 


3 


00 


-0 


39 


24 


29 


4 


21 


Fe I 


6269.55 





86 


-1 


33 


70 


86 


15 


98 


Fe I 


5283.63 


3 


24 


-0 


52 


17 


20 


3 


20 


Fe I 


5324.19 


3 


21 


-0 


10 


25 


30 






Fe I 


5328.05 





92 


-1 


47 


68 


81 


13 


82 


Fe I 


5371.60 





96 


-1 


64 


56 


66 


12 


77 


Fe I 


5383.38 


4 


31 


+0 


65 


17 


20 






Fe I 


5397.14 





92 


-1 


98 


47 


66 


5 


32 


Fe I 


6406.78 





99 


-1 


86 


44 


62 


9 


61 
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Table 3 — Continued 



S i 
pecies 


\ 




X 


1 f 
iog 9/ 


EWp o 


EW-p t 


EWs o 


EWs t 




A 


eV 








Fe I 


6429.71 





96 


-1.88 


46 


66 


10 


59 


Fe I 


5434.53 


1 


01 


-2.13 


36 


44 


7 


40 


Fe I 


5446.92 





99 


-1.91 


39 


49 


20 


106 


Fe I 


5501.48 





96 


-3.05 


10 


12 






Fe I 


5506.79 





99 


-2.79 


15 


18 


5 


29 


Fe I 


6669.63 


3 


42 


-0.60 


11 


14 


3 


17 


Fe I 


6672.86 


3 


40 


-0.28 


16 


19 






Fe I 


6686.77 


3 


37 


-0.10 


21 


26 


5 


30 


Fe I 


6616.66 


3 


33 


+0.05 


24 


28 


5 


27 


Fe I 


6136.62 


2 


45 


-1.41 


16 


19 






Fe I 


6191.57 


2 


43 


-1.42 


13 


16 






Fe I 


6230.74 


2 


56 


-1.28 


16 


19 


3 


17 


Fe I 


6393.61 


2 


43 


-1.68 


13 


16 






Fe I 


6400.01 


3 


60 


-0.29 


13 


16 






Fe I 


6430.86 


2 


18 


-1.96 


10 


12 






Fe I 


6494.98 


2 


40 


-1.24 


22 


26 


12 


66 


Fe II 


4178.86 


2 


57 


-2.48 


20 


23 






Fe II 


4233.17 


2 


58 


-2.00 


37 


44 


3 


19 


Fe II 


4508.28 


2 


84 


-2.22 


16 


19 






Fe II 


4666.89 


2 


82 


-2.28 


12 


14 


5 


33 


Fe II 


4923.92 


2 


89 


-1.32 


48 


56 


8 


47 


Fc II 


5234.62 


3 


22 


-2.05 


17 


20 






Fc II 


5275.99 


3 


20 


-1.94 


15 


18 






Co I 


3502.29 





43 


+0.07 


52 


61 


11 


81 


Co I 


3506.33 





51 


-0.04 


39 


46 


9 


65 


Ni I 


3472.56 





11 


-0.82 


55 


64 


11 


79 


Ni I 


3483.78 





27 


-1.12 


60 


67 






Ni I 


3500.86 





17 


-1.37 


54 


62 


12 


89 


Ni I 


3510.33 





21 


-0.65 


73 


84 


12 


90 


Ni I 


3515.07 





11 


-0.26 


70 


81 


17 


124 


Ni I 


3519.76 





27 


-1.42 


36 


41 


5 


37 


Ni I 


3524.54 





03 


-0.03 


84 


97 


16 


118 


Ni I 


3566.39 





42 


-0.27 


64 


74 


12 


84 


Ni I 


3697.71 





21 


-1.09 


44 


61 


12 


87 


Ni I 


3619.40 





42 


-0.04 


82 


96 


7 


50 


Sr 11" 


4077.71 





00 


+0.17 


(86) 


(101) 


(16) 


(110) 


Sr 11" 


4215.52 





00 


-0.17 


(80) 


(94) 


(16) 


(99) 


Y II 


3611.04 





13 


+0.01 










Y II 


3710.29 





18 


+0.46 










Y II 


3774.33 





13 


+0.21 










Zr II 


4208.99 





71 


-0.46 










Ba II* 


4564.03 





00 


+0.17 


(92) 


(108) 


(13) 


(82) 


Ba II" 


4934.10 





00 


-0.16 


(82) 


(98) 


(10) 


(62) 


Ba II 


5853.69 





60 


-0.91 










Ba II" 


6141.73 





70 


-0.08 


(58) 


(71) 


(8) 


(46) 


Ba II" 


6496.91 





60 


-0.38 


(63) 


(66) 


m 


(36) 


La II 


3988.51 





17 


+0.21 










La II 


3996.74 





17 


-0.06 










La II 


4086.71 





00 


-0.07 










La II 


4123.22 





32 


+0.13 










Co II 


4137.65 





52 


+0.44 










Nd II 


4021.33 





32 


-0.10 










Nd II 


4061.08 





47 


+0.55 










Nd II 


4462.98 





66 


+0.04 










Eu II 


3819.67 





00 


+0.51 










Eu II 


3907.11 





21 


+0.17 










Eu II 


3971.97 





21 


+0.27 










Eu II 


4129.72 





00 


+0.22 










Eu II 


4205.05 





00 


+0.21 










Yb II 


3694.17 





00 


-0.30 










Pb I 


4057.81 


1 


32 


-0.17 











^ EW values are given for illustration only, and were not used in final abundance 
determinations. 
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Table 4. Individual Abundance Results for Primary and Secondary Stars 



Species 


log Cq 


log Cp 


ap 




log Cs 






[X/B.]p 




[X/Fe]p 


[X/Fe], 


Method 


CH 


8.70 


7.64 


0.10 




7.44 


0.40 




-1.06 


-0.86 


+ 1.35 


+1.15 


syn 


Fe I 


7.52 


5.13 


0.16 


104 


5.11 


0.44 


83 


-2.39 


-2.41 


+0.02 


0.00 


EW 


Fe II 


7.52 


5.09 


0.12 


7 


5.11 


0.44 


3 


-2.43 


-2.41 


-0.02 


0.00 


EW 


Nal 


6.33 








3.9 


0.3 


2 




-2.4 




0.0 


syn 


Mgl 


7.58 


5.53 


0.20 


6 


5.52 


0.15 


6 


-2.05 


-2.06 


+0.36 


+0.35 


EW 


All 


6.57 


3.26 




1 


3.28 




1 


-3.31 


-3.29 


-0.90 


-0.88 


EW 


Cal 


6.36 


4.38 


0.13 


8 


4.03 


0.29 


6 


-1.98 


-2.33 


+0.43 


+0.08 


EW 


Ti II 


4.99 


3.19 


0.26 


9 


3.14 


0.75 


4 


-1.80 


-1.85 


+0.61 


+0.56 


EW 


Cr I 


5.64 


3.36 


0.30 


8 


2.81 


0.53 


6 


-2.28 


-2.83 


+0.13 


-0.42 


EW 


Cr II 


5.64 


3.45 


0.10 


2 


3.64 




1 


-2.19 


-2.00 


+0.22 


+0.41 


EW 


Ni I 


6.25 


3.92 


0.24 


10 


3.86 


0.62 


9 


-2.33 


-2.39 


+0.08 


+0.02 


EW 


Sr II 


2.90 


1.05 


0.09 


2 


0.93 


0.29 


2 


-1.85 


-1.97 


+0.56 


+0.44 


syn 


Ba II 


2.13 


1.17 


0.08 


1 


1.09 


0.28 


1 


-0.96 


-1.04 


+1.45 


+1.37 


syn 
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Table 5. Abundance for the CS 22964-161 System 



Species 


log e© 


log e 


cr 


# 


[X/H] 


[X/Fe] 


Method^") 


"Pq t 
re i 


i .oZ, 


-TO. lO 


n 91 

U.Zi 


iU4 


9 'iQ 
— Z.Oil 


_Ln 09 

+U.UZ 


est 


T?Q TT 

re ii 




no 
+O.Uy 


n 1 7 

U.l ( 


( 


9 /I Q 


n n9 
— U.UZ 


est 


T i T 




-|-z.uy 


n 9n 


1 
i 






est 




o. ( U 


_l7 c;n 
+ ( .ou 


u.io 


1 
i 


1 9n 

— i.ZU 


_Ll 91 
-t-i.Zi 


syn 


1\Tq t 
iNa i 


u.oo 


-L'? Q 

-i-o.y 


u.o 


A 
4 


9 A 
— Z.4 


n n 
u.u 


syn 


A/Trr T 


7 c^Q 
( .Oo 


"t-O.OO 


n 9n 





9 nc; 
— z.uo 


_i_n 'ifi 


syn 


Al T 


0.0 ( 




n 9n 


1 
i 


— O.Oi 


n an 
— u.yu 


est 




( .00 


J-K 98 
-t-O.Zo 


u.zo 


1 
i 


9 97 
— Z.Z 1 


-LD 1 /I 
-|-U. i4 


syn 




o.oo 


^^9 


U. iO 


8 

o 


9 D/l 


-LD ^^7 
-|-U.O ( 


est 


Ti TT 
1 1 ii 


/I QQ 


-|-0. iO 


U.04 


Q 

y 


1 81 
— i.oi 


-ld fin 


est 


T 

or i 


0.04 


_lQ 97 


U.o4 


Q 
O 


9 Q7 
— Z.o 1 


+U.U4 


est 


Pv TT 

or ii 


0.04 


_lQ /18 


n 1 9 

U. iZ 


Z 


9 1 
— Z. iO 


_Ln 9c; 

+U.ZO 


est 


A/Tn T 
iVill i 


o.oy 


_l9 97 


U.iO 


Q 
O 


1 9 
— O. iZ 


n 71 

— U. ( i 


syn 


iNl i 


O.ZO 


+o.yl 


U.oU 


1 n 
iU 


9 Q/1 

— Z.o4 


_Ln n7 

+U.U ( 


est 


Sr II 


2.90 


+1.01 


0.11 


2 


-1.89 


+0.52 


syn 


Y II 


2.24 


+0.12 


0.10 


3 


-2.12 


+0.29 


syn 


Zr II 


2.60 


+0.88 


0.25 


1 


-1.72 


+0.69 


syn 


Ba II 


2.13 


+0.98 


0.16 


5 


-1.15 


+1.26 


syn 


La II 


1.13 


-0.21 


0.10 


4 


-1.34 


+1.07 


syn 


Ce II 


1.55 


+0.43 


0.25 


1 


-1.12 


+1.29 


syn 


Nd II 


1.45 


+0.05 


0.12 


3 


-1.40 


+1.01 


syn 


Eu II 


0.52 


-1.20 


0.26 


5 


-1.72 


+0.69 


syn 


Yb II 


1.08 


-0.69 


0.25 


1 


-1.77 


+0.64 


syn 


Pb I 


1.85 


+1.63 


0.20 


1 


-0.22 


+2.19 


syn 



Abundance method: "est" means that the system abundance 
was estimated from the individual primary and secondary abun- 
dances given in Table HI "syn" means that the system abundance 
was computed with synthetic spectra of features in the syzygy spec- 
trum. 
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Table 6. Carbon-Rich Star Periods and Velocity Ranges 



Group 


P^^i'ange 


Median P 


K'^^i'ange 


Median K 




days 


days 


km s~^ 


km s~^ 


Barium-CH giant 


80 - 4390 


1352 


2.5 - 12.0 


6.0 


CH subgiant 


880 - 4140 


1930 


3.3 - 7.2 


4.6 



orbital period 
(^^velocity semi- amplitude 
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Table 7. Predicted mass fraction and molecular weight 



Mass Fraction M^l.SM^ M=1.5M 



At last Third Dredge-Up 



X 


0.675 0.535 


Y 


0.301 0.403 


Z 


0.0238 0.062 


1^ 


0.630 0.713 


Mass average of the winds 


X 


0.700 0.585 


Y 


0.284 0.367 


Z 


0.0155 0.0483 


1^ 


0.616 0.781 



